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M e thod for d e t e rmining th e humidity and d e nsity of a di e l e ctric mat e rial 
METHOD FOR DETERMINING THE HUMIDITY AND DENSITY OF A 
DIELECTRIC MATERIAL 



Background of the Invention 

The invention relates to a method for determining the humidity and/or density 
of a dielectric material in a resonator filled with the material, comprising a 
sender and a receiver. 



The dielectric properties of a material, described by the complex relative 
dielectric constant £ r = s[ - ie" 9 can be affected by humidity and density in case 
of porous materials. Humidity and density thus modify the scalar parameters 
resonant frequency fir and resonator quality Q of a material-filled resonator in 

contrast to those of the empty, air-filled resonator (fr 0 an d Qo) in such a way 
that 

f — 1 1 fi r 

rm ~ and Q^ = Q^ + T, wherein 1 ^ £ r £ eL» and 0 ^ £ r < C*, wherein 

£ rmax and £ r ma x are the maximum values resulting for the respective material 
from the assigned humidity range and density range. 

In prior art, different methods are known for measuring the humidity or density 
of granular materials, in which the resonance behavior of a matter-filled 
resonator is used. 



For example, from US 566 60 61 Patent No. 5,666,061 a method is known for 
measuring the humidity in granular materials by means of microwaves. There, 
one electronically reacts to a threshold while bidirectionally sweeping a 
frequency range. From the chronological properties of pulses created in such a 
way the components of the dielectric constant are deduced. 
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All methods are relatively slow and have either a dependency of the density of 
the material or have relatively large errors in the determination of the humidity. 

The invention is underlied by the problem to specify a method of the type 
initially mentioned, by which a fast, accurate and density-independent 
determination of the humidity is possible. 

Summary of the Invention 

According to th e inv e ntion, th e probl e m is oolv e d by a m e thod comprising tho 
attribut e s giv e n in claim 1 . 

Advantag e ous e mbodim e nts ar e giv e n in th e d e p e nd e nt claims. 

For the purpose of the invention, any arbitrary quantity describing the 
resonance width can be considered as a bandwith in the following, wherein 
appropriate adjustments have to be provided in the respective definitions, in 
particular of the threshold values, and equations. 

By digitally recording the resonance curve A a fast acquisition of the measuring 
values is possible. Thus A a close chronological-spatial assignment of the 
humidity determined from these values to material dynamically guided through 
the resonator is possible. For this purpose, the invention provides for that the 
sender emits a signal; a resonance curve of the filled resonator is sweeped, 
wherein respective related signal strength values of the receiver signal are 
measured at different frequencies; the resonant frequency and the bandwidth 
are determined for the filled resonator from the measured points; and the 
humidity and/or density of the material is calculated by solving a second 
system of equations comprising the resonant frequencies and the bandwidths of 
the empty and of the filled resonator and known calibration coefficients of the 
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resonator. This method enables a density-independent determination of the 
humidity of the material. Besides, the density of the material can be acquired 
with little effort. 

A preferred embodiment provides for that, from the points for determining the 
bandwith of the filled resonator, either the quantities resonant frequency, 
resonator quality and resonance maximum are determined and the bandwidth is 
calculated therefrom, or cut-off frequencies are determined and the resonant 
frequency and the bandwidth are calculated therefrom. 

In another embodiment a lower threshold value is calculated and a second 
sweeping pass with smaller step sizes is performed in that range in which the 
signal strength values are higher than the threshold value. By a two-pass 
procedure the accuracy can be significantly increased and, nevertheless, the 
required time can be kept short, in particular if, during the second sweeping 
pass, a decreased step size is used in a range of the resonance peak only. 

Advantageously, sweeping the resonance curve is performed in equally spaced 
steps. The simplest and fastest form of sweeping consists of equally spaced 
steps. Variable step sizes in the second sweeping pass can shorten its duration. 

Preferably, the sender is operated using a constant strength. The measuring 
values at the receiver can be used without adjusting or scaling them in case of a 
constant signal strength. 

In a possible embodiment^ the cut-off frequencies of the resonator are 
determined by determining the point having the highest receiver signal strength 
value, and, starting from this point, calculating a threshold value; determining 
two respective proximate points for positive and negative slope sections, the 
signal strength values of these points lying below and above the threshold 
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value; calculating first and second cut-off frequencies therefrom by 
respectively interpolating between the proximate points. Determining the cut- 
off frequencies by interpolating between point pairs surrounding a threshold 
value as an initial parameter for determining the humidity and/or density is a 
fast and simple method. 

An advantageous embodiment thereby provides for that the threshold value 
corresponds to an attenuation of 3 dB in relation to the highest signal value. If 
the threshold value is chosen corresponding to an attenuation of 3 dB, starting 
from the maximum signal strength value, the equations to be solved obtain a 
very simple form. 

In an alternative embodiment the quantities resonant frequency, resonator 
quality and resonance maximum of the resonator are determined by arbitrarily 
and/or randomly selecting three points and solving a first system of equations 
for these quantities, the system consisting of three equations of an analytic 
resonance curve that are valid for the three points. The resonant frequency, the 
resonator quality and the resonance ampitude can be directly determined as the 
initial parameters for determining the humidity and/or density even faster and 
with less error than by interpolation by solving a first system of equations 
consisting of three equations and thus being completely determined, wherein 
for each point out of a group of three one of the three equations is valid, the 
three points being selected from the present set. In particular, noisy resonance 
curves can be analyzed more exactly this way. 

In another alternative embodiment the quantities resonant frequency, resonator 
quality, resonance maximum of the resonator are determined by arbitrarily 
and/or randomly selecting a set of points whose number is an integer multiple 
of three and at least six, and splitting up the point set into three equally sized 
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groups; for each combination of three points, wherein each point comes from a 
different group, solving a first system of equations for these quantities, the 
system consisting of three equations of the analytic resonance curve valid for 
these three points; and creating the average for each quantity from the values 
calculated at the combinations. Even more exact values than with three points 
can be obtained by creating several triple-groups from the present point set and 
averaging the initial parameters obtained thereby over all groups. 

In an advantageous embodiment, as a condition for arbitrarily and/or randomly 
selecting the points, the signal value of a point to be selected is higher than the 
highest signal value attenuated by 3 dB. In both procedures using three points 
or a multiple thereof, preferably only those points from the resonance curve are 
selected whose signal strength is higher than the maximum value of all 
measurement values, attenuated by 3 dB, as this way only significant values are 
used. 

The system of equations to be solved is preferably chosen in such a way that it 
describes, in a good approximation, the correlation of humidity and density 
with a variation of resonant frequency and resonator quality or, respectively, 
with a variation of resonant frequency and bandwidth. 

Advantageously the second system of equations is non-linear. 

In a preferred embodiment the sweeping by means of the sender is performed 
up to the microwave area. 

Advantageously, voltage values or current values of the receiver are used for 
measuring the receiver signal Preferably, the electrical voltage rising at the 
receiver serves for acquiring the resonance curve, because it is measurable 
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easily and without back coupling. However, the current in a receiver circuit can 
be measured, too. 

In the following, the denomination Ui i s used for the measured quantities. 
However, in doing so, it shall not be implied to use the voltage only. 

The resonance curve to be reconstructed from the discrete measuring points is 
an approximation as-good-as-possible to the real resonance curve only. It is 
characterized by the course of the signal amplitude U a t the receiver against the 

supplied frequency f, wherein U r i s the resonance maximum and Q is the 
quality of the resonance: 

The cut-off frequencies fa and f b are those frequencies at which the signal 
strength exceeds and falls below a defined value, respectively. For this purpose, 
preferably the value corresponding to an attenuation of the maximum by 3 dB 

is chosen: a " h ~~ Ji . Between the cut-off frequencies and the curve 
parameters resonant frequency f r and resonator quality Q there are the 

correlations: fr = V f a • fi> , Q = f b _f a , The distance between the cut-off 
frequencies is defined as the bandwidth = fb - f a = "q 



In the following, the invention is described in further detail using examples of 
embodiments. 

Brief Description of the Drawings 
For this purpoGo, 

Figure 1 shows a sweeped resonance curve[[,]]; 
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Figure 2 shows a resonance curve sweeped in a two-pass procedure[[,]]; 
Figure 3 shows a first way for determining the resonance parameters[[,]]; 
Figure 4 shows a second way for determining the resonance parameters[[,]]; 
Figure 5 shows a third way for determining the resonance parameters; 
5 and 

Figure 6 shows a two-dimensional representation of calibration curves. 
Detailed Description of the Invention 

In Figure 1 , two resonance curves are depicted above each other like they have 
10 been recorded using the method according to the invention; the right resonance 

curve at an empty, the left resonance curve at a material-filled resonator. 

By sweeping (wobbling) the resonator between the start frequency fstarti and 
the stop frequency fstopi 4 the resonance curve is acquired in discrete steps. The 
15 start frequency is calculated from the maximally shifted frequency 

V c "»« , the maximally changed quality Qm ~ Qo + b ^ and the normalized 

n - Ua <- -t- f - - frm I±=SL , Tt fr ™ O 1-a 2 " cO 

voltage ratio u — : IstartI ~ 2Q m V a> ^^2Q m ) • a 2 + ir m The 
stop frequency, using fr = fr 0 , results as: 

20 

The sweeping speed depends, among others, on the number n i of sweeping 
points and can be increased by a two-pass procedure as depicted in Figure 2. 
For this purpose, first a sweeping is performed using a smaller number of 

sweeping points and, hence, a larger frequency step size Atl ~ n,-i and, in 
25 a second iteration, after previously determining the start and stop frequencies 

fstarfl =f\u*U a A/>0, fstopl =f\(J*U a A/<0 or fstard ~f\ U*U a A U<U ry 
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fstopi-f\u*u a t\U>Ur, another sweeping is subsequently performed using a 

* ~ fstop2~fstart2 

smaller frequency step size At2 " n 2 -i between these frequencies . 

From the sweeped, measured resonance curve the resonator parameters 
f r , Q and U r can b e determined, for example, according to the fafb procedure, 
the three-points procedure or the 3k-points procedure. Other procedures are 
possible, too. 

Figure 3 shows the fafb procedure. It is based on directly determining a first 
and a second cut-off frequency fa, fb from the measured resonance curve. For 
this purpose, first the point (also called sweeping point) 
(f = fmax/U = Umax)having the highest voltage U=U max [ s determined for 

j j Umax 

calculating a 3 dB threshold line 3dB " JT . Subsequently, the respective two 
points in the immediate proximity of the 3 dB threshold line are determined. By 
linearly interpolating between the sweeping points al and a2 as well as bl and 
b2 the first and the second 3 dB cut-off frequency fa, fb are obtained: 

Umax 1 1 Umax » , 

jY -Uai -Jf — U b | 

fa = fal + Ua2-U a , ( f a2 ~ fa) ) ? fb = fbl + U b2 -U b i ( f b2 ~ fbl ). 

From the equations fr = V f a • f b , Q = "^fT ? f r and Q can be calculated. The 
resonance voltage U r then results as: 

Calculating the cut-off frequencies is more or less error-prone because of the 
discretization of the resonance curve - on one hand, by the determination of the 
3 dB threshold line (U a = U max /72 < U r /V2 ) 

derived from U max and, on the 
other hand, by the linear interpolation between the sweeping points al, a2, bl, 
b2. A higher sweeping rate indeed reduces the errors, but also increases the 
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sweeping time and thus reduces the sweeping speed. At best, a non-linear 
interpolation is possible in this procedure in order to improve accuracy. 

In the three-points procedure shown in Figure 4, three sweeping points 
U,,U 2 ,U 3 

are selected arbitrarily or by a random generator preferably above 
the line U s = s • Umax with S < 7T. 



By solving the first system of equations Gl : 

u ' = MU) 2 , U2= MHO 2 , U3= 

the resonator parameters/', Q an d Ur are obtained: 

_ / u?(ul-ua)f?-rtJi(uf-ui)fi-rtjj(ul-uf)fi 

r 4 uf(uj-uj) ufrufofo u]qj-u^) 
1 f i + f i + 1 



o= a^j 



The advantage of the three-points procedure in relation to the fafb procedure 
lies in that no errors occur because of the discretization of the resonance curve, 
and that only few points are needed for sweeping the resonance curve A and thus 
a high sweeping speed is achieved. 

If the resonance curve is noisy, error-prone resonator parameters are 
determined by the three-points procedure, too. In order to reduce the impact of 
the noisei more than three sweeping points are selected, namely a 
number^-k, k = 2,3,4, ...^ anc j m sp ij t up i nt0 three groups Ml, M2, M3, 
having k points per group above the line U s = s • t Figure 5 shows this 3k- 
points procedure. For all k 3 possibilities of combining one respective point 
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from each of the three groups, the resonator parameters are determined from 
the first system of equations Gl according to the three-points procedure 

described above, and the k 3 partial results are averaged subsequently: 

fr = f *S K Q = f zft Qj, u r = tJ- if, u rj . 

In order to determine humidity and density from determined resonator 
parameter the used arrangement is calibrated using material of known 
humidity and density. The calibration is advantageously performed out of the 
regular operation. 

In the calibration, the electrical resonator values resonant frequency fr, 
resonator quality Q and bandwidth BW(BW = f r /Q) 

are respectively assigned 

to the material quantities humidity W and density Q, This assignment between 

the material values and the electrical values is non-linear and can be described 

for a defined range of humidity and density with sufficient accuracy, for 

example, by the following second system of equations G2: 

Afr = HiQ 1 ¥ 2 + Hi QV 2 + H2Q 2 ¥ + bfn Q¥ + Cf«Q 2 + Cf rI Q 

ABW = a B w 2 £>V 2 +a B w,ey 2 + b B w 2 eV +bBw,ey + cbw 2 £ 2 +c B w,e 

Afr and ABW are the differences of resonant frequency and bandwidth, 
respectively, between the empty and the material-filled resonator: 
Af r = f ro - f rmj ABW = BW m - BWq 

Now, it is an object of the calibration to determine the twelve calibration 
coefficients a fn> a fr2> Km bf r2 , Cf rl , Cf r2 , a B w,, aBw 2 , bBw,, bBw 2 , Cbw,, cbw 2 
for the respective material from a sufficient large number of calibration values 
(Af r and ABW with related material values V und Q). 
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For this purpose, the measuring values of the resonant frequency difference Af r 
and the bandwidth difference ABW are assigned to the humidity and density 
values that have been determined by appropriate reference methods. These 
calibration values characterizing the resonant frequency and bandwidth 
depending on humidity and density are the basis for determining the calibration 
coefficients. The regression of the calibration values of equal humidity is 
performed in the form of Af r = d fr2 ^ 2 + d fr ^ a nd Af r = d B w 2 £> 2 + d B w,e, 
wherein the regression curves have to run through the origin, as Af r and ABW 
are zero for the empty resonator (air having Q = 0). The regression provides 
calibration points of equal density depending on the humidity. Using these, 
another regression of the form Af r = af r ^ 2 +bf r ^ + Cf r a nd 
ABW= a B w^ 2 + b B w^ + c B w is performed. The values determined this way 
for the regression coefficients a fr> bf r » Cf r an d a B w, b B w, c BW are plotted 
against the density and, therefrom, the calibration coefficients are determined 
by a quadratic regression. 

Figuf -Figure 6 shows the results of the calibratrion in the form of a set of 
curves which also serves for analysis. 

From the measured resonant frequency difference values Af r and bandwidth 
difference values ABW humidity content ¥ and density Q are calculated for the 
respective material by solving the above second system of equations G2. 

In doing so, two real and two imaginary roots [[do]] result. From the course of 
the calibration curves in the Af r _ABW diagram it can be determined if there is 
only one real solution in the humidity and density range of interest. In the Af r . 
ABW diagram the bandwidth difference is represented in dependency of the 
resonant frequency difference for curves of equal density and humidity. If the 
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course of these curves, characterized by the points A, B, C und D in Figure 6, is 
continuous and unique in the humidity and density range of interest then only 
one real solution exists in this range. 

For solving the second system of equations G2 an iterative procedure is 
appropriate. For this purpose, the second system of equations G2 is solved for 
¥: 

_ _ b fr2 g+b fr i f? b fr2 g+b fr i \ 2 ~ Cf r2 g 2 +Cf r i-Af ro 
t 2a fr2 <H-2a fr i£ + ^ ^2a fr 2^+2a frl e J a fr2 e 2 +a fr ie 

_ b BW 2 (?+bBW 1 \f bBW 2 g+bew 1 \ 2 cbw 2 ^ 2+ cbw,-ABWo 
V — 2aB\v 2 ^+2aBW 1 p y ^2aBW 2 CH-2aBW,e J _ aBW 2 0 2 +a B w,£> 

From the intersection of both equations in a humidity-density diagram the 
sought values for ¥ and Q are obtained. 
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List of reference numerals 

fr Resonant frequency in the general case 

^0 Resonant frequency of the empty resonator 

5 fr- Resonant frequency of the filled resonator 

fstarti , fstopi Start and stop frequencies for first sweeping pass 

Umax Highest signal strength value measured 

a, s Threshold value factors 

fstart2, fstop2 Start and stop frequencies for second sweeping pass 

10 fmax Frequency at which the highest signal strength 

value is present 

(fai/Uai), (fa2/Ua2) Proximate points of the first cut-off frequency 

(fbi/Ubi), (fb2/Ub2> Proximate points of the second cut-off frequency 

(fi/Ui) ... (f3/U 3 ) Three selected points 

15 M 1 , M2, M3 point groups 

(fii/Un) ... (f i 4 AJ 14) Elements of point group Ml 

(f2i/U 2 i) ... (f24/U 2 4> Elements of point group M2 

(f 3 i/U 3 i) ... (f34/U34> Elements of point group M3 

A,B>C,D Limits of the humidty and density range of interest 
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